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I. INTRODUCTION 



The quantum and semiclassical theory of Dirac particles in the gravitational field was 
studied in the numerous papers (see Refs. jl| to mention but a few). The main attention 
was usually paid to the dynamics of fermions in the weak gravitational fields when the com- 
ponents of the spacetime metric tensor do not significantly deviate from the flat Minkowski 
metric. The analysis of the strong gravitational fields is of interest for the investigation 
of the processes near black holes and massive compact astrophysical objects. They may be 
also applied to describe the spaces of variable geometry and dimension discussed recently jj- 
4|. Yet another application corresponds to the fermions (quarks) motion in rotating media 
which may be realized j^J in heavy-ion collisions. 

In this paper, we apply the methods of the Foldy-Wouthuysen (FW) transformation pre- 
viously used for the analysis of the spin dynamics in weak static and stationary gravitational 
fields for the case of arbitrary strong gravitational fields. We use the isotropic spa- 
tial coordinates and systematically investigate, in particular, the spin dynamics of Dirac 
fermions. Specifically, we derive the general result for the angular velocity of spin precession 
and compare this quantum mechanical result with classical predictions. 

The paper is organized as follows. In Sec. [TT] we derive the general Hermitian Dirac 
Hamiltonian for a charged particle in an arbitrary curved spacetime interacting with the 
electromagnetic fields. In this general framework, we further specialize to the dynamics of 
the spin and momentum of the particle in a spatially isotropic metric. Sec. II III presents the 
FW Hamiltonian and operator equations of motion. The equations describing the dynamics 
of a classical spin are obtained in Sec. IIVI and we compare them with the corresponding 
quantum equations, testing the equivalence principle. The validity of the latter as well as 
the similarity of the spin interactions with the gravitational and with the electromagnetic 
fields is also manifested in the analysis of squared covariant Dirac equation performed in 
Sec. IVI1 The evolution of the helicity in the strong fields is discussed in Sec. [V] The results 
obtained are summarized in Sec. IVHI 

We denote world indices by Latin letters i,j,k,... = 0,1,2,3 and reserve first Greek 
letters for tetrad indices, a, (3, . . . = 0, 1,2,3. Spatial indices are denoted by Latin letters 
from the beginning of the alphabet, a,b,c,... = 1,2,3. The separate tetrad indices are 
distinguished by hats. As usual, A and * denote the exterior product and the Hodge operator, 
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respectively. 



II. SPINNING PARTICLE IN CURVED SPACETIME 
A. General spacetime metric 

A classical spinning particle is characterized by its position in spacetime, x 1 {t) which 
is a function of the proper time r, and by the 4- vector of spin S a . The 4- velocity of 
a particle U a = efdx l /dr is normalized by the condition g a pU a U 13 = c 2 where g a p = 
diag(c 2 , —1, —1, —1) is the flat Minkowski metric. To describe spinning particles in flat and 
curved spacetimes (as well as in arbitrary curvilinear coordinates), we use the tetrad e°. 
When the gravitational field is absent, one can choose the Cartesian coordinates and apply 
the coincidence of the holonomic orthonormal frame with the natural one (e" = 5f). For an 
arbitrary spacetime metric, g a pe^e^ = gij. Let t be the time and x a (a = 1,2,3) be spatial 
local coordinates. The general form of the line element of an arbitrary gravitational field 
can be given by 

ds 2 = V 2 c 2 dt 2 - 5^W^ c W l d (dx c - K c cdt) (dx d - K d cdt). (2.1) 

The functions V and K a , as well as the components of the 3x3 matrix W a b may depend 
arbitrarily on t, x a . The total number 1 + 3 + 9 = 13 of these functions is larger than the 
number of components of the spacetime metric. However, it is obvious that the structure of 
(12. ip allows for a redefinition W a b — > L a zW C b with the help of an arbitrary local rotation 
L a ^{t, x) G SO(3). Taking this freedom into account, we end up with exactly 10 independent 
variables that describe the general spacetime metric. 

The off-diagonal metric components go a are related to the effects of rotation, as is well 
known. They arise when the functions K a are nontrivial. For example, the exact metric of 
the flat spacetime seen by an accelerating and rotating observer is obtained for the special 
case [lo| 

V = l + ^f, W" b = 5 b 1 , K a = --(uxr) a , (2.2) 

where a describes acceleration of the observer and oj is an angular velocity of a noninertial 
reference system. Both are independent of the spatial coordinates, but may depend arbitrar- 
ily on time t. The slowly rotating massive body produces the gravitational field of similar 
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form; far from the source the Lense-Thirring ll| metric is described by 

V = V(x), W\ = 51 W(x), K a = ^-e abc u b {x)x c . (2.3) 

The functions V(x a ), W(x a ), uj(x a ) can be recovered from the Kerr metric in the isotropic 
coordinates in the limit of r — > oo: 

V =( 1 ~B + + 0( a ^), (2.4) 

^(ll^'l ^-^ ^M (2.5) 
and if a is given by the same equation ( 12.31) with 

u= 2 J^ a ( 1 - 3 ± + 2 24) +0 (a 3 r-% (2.6) 



V ,y /J^y 2 

Here r := \Jx ■ x and n = r/r. The constant vector a = (0, 0, a) is the rotation parameter 
of the Kerr solution and a ■ n = az/r. Also, fi = GM/c 2 ; the total mass M and the total 
angular momentum J = Mca define the Kerr black hole uniquely. These equations are 
obtained from the Arnowitt-Deser-Misner form [3] of the Kerr solution performed earlier 



by Hergt and Schafer 13[ after dropping the terms violating the isotropy. 



In the weak field approximation, we have studied the dynamics of quantum and classical 



spin for the cases d J and $ 3]) in Ref. (o | - 



B. Dirac fermions 



In order to discuss the Dirac spinors, we need the orthonormal frames. The preferable 



choice [9] is the Schwinger gauge: 

ef = V 5°, ef = W\ {5\ - cK b 6?), a = 1, 2, 3. (2.7) 

Tetrad (12. 7p is characterized by the condition e° = 0,a = 1,2,3. The same is automatically 
true for the inverse tetrad: 

e\ = 1 (<f + 5\cK a ) , el = 5\W b ^ a = 1,2, 3, (2.8) 

where the inverse 3x3 matrix, W a cW C b = $£, is introduced. 

The covariant Dirac equation for spin-1/2 particles in gravitational and electromagnetic 
fields has the form 

{ihj a D a - mc)^ = 0, a = 0,1, 2, 3. (2.9) 



The Dirac matrices 7° are defined in local Lorentz (tetrad) frames. The spinor covariant 
derivatives are given by 

D a = e\D h Di = di + ^Ai+ % -a afi Y iafi . (2.10) 

Here = — r/ a are the Lorentz connection coefficients, a alB = - (7°7 /3 — 7 /3 7°)- The 
Dirac particle is characterised by the electric charge q, and A\ is the 4-potential of the 
electromagnetic field. 

Eqs. ( 12.91) . (12.101) show that the gravitational and inertial effects are encoded in coframe 
and connection (see the relevant discussion in Refs. [14J, [15[ and references therein). For the 
general metric (12. ip with the tetrad (12. 7p we find explicitly 

r, S o = C ^ b ad h Ve?-^Q^el (2.11) 
r i{ £ = ^ Qm e? + + + C^) ef. (2.12) 



Here 



= gaeW% ( -W^ d + K e d e W^ d + W\d d K e ) , (2.13) 



c 

C af = W d ^ b d [d W% = -G£, C afe = g^C z {. (2.14) 

The dot ' denotes the derivative with respect to the time t. As one notices, C^ c is nothing 
but the anholonomity object for the spatial triad W a b- The indices (that all run from 1 
to 3) are raised and lowered with the help of the spatial part of the flat Minkowski metric, 
9ab = -Sab = diag(-l, -1,-1). 

Dirac equation can be derived from the action 

d*x£, C = y/^jL (2.15) 

with the Lagrangian 

L = — (^ a D a ty - DjJ^) - mctftf. (2.16) 

As it is well known, the naive Hamiltonian form the Schrodinger form of the Dirac equation 
is not Hermitian. The most straightforward way to derive a Hermitian Hamiltonian is to 
redefine the wave function. Substituting the tetrad and connection into the Lagrangian 
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density, we find explicitly 



dtd x 
Hi 



i h 



2 



- 4 v^^r/^tS^ + -y/^ge^T^^ 



(2.17) 



A direct check shows that the Schrodinger equation derived from this action has a non- 
Hermitian Hamiltonian. This problem disappears if we define a new wave function by 



(2.18) 



Such a non-unitary transformation may be also interpreted in the framework of the pseudo- 
Hermitian quantum mechanics [16|, ll7j (cf. also [18]). Substituting (12.181) into (I2.17p . and 
recalling ( 12. 7ft and ( 12. 8ft . we find the action 



ih 



— Wa a d b ^ - d b ^a a ij) -qA h ^a a i) 

he 

--(S a ^E^ + T^W) 



(2.19) 



Here V = e° , T h a 
T = 



v/=?4 = and 



^ afec 



be 



(2.20) 



abc abc" 1 

For the static and stationary rotating configurations, the pseudoscalar invariant vanishes 

'ab? 



^e a6c C s j 3 _^= 0), and thus the corresponding term was absent in the special cases considered 



earlier 



9[. But in general this term contributes to the Dirac Hamiltonian. 



Variation of the action (12.19P with respect to the rescaled wave function yields the Dirac 

dip 

equation in Schrodinger form ih^— = Hib. The corresponding Hermitian Hamiltonian reads 

at 



-H = (3mc 2 V + g$ + - (n b T l \a a + a a F b a ii b ) 



'K-ir + ir-K) 



he 



(T 75 + H-S). 



(2.21) 



q q 

The kinetic momentum operator 7Tj = ihdi A^ = Pi accounts of the interaction 

c c 

with the electromagnetic field Ai = (<3>, A a ). 

Hamiltonian f)2.2ip is one of our central results and covers the general case of a spin-1/2 
particle in an arbitrary curved spacetime. Remarkably, the form of the Hamiltonian remains 
the same even when the connection is non- Riemannian. If the torsion T a = dfi a + V p a A $ p 
is nontrivial, we just have to replace 

~ a ^~ a + T a , (2.22) 

where T a d a = *(i9 a AT a ) is the axial covector of torsion (i.e., the totally antisymmetric piece 
of torsion). It is worthwhile to mention that the recent discussion Ij]] of the Dirac fermions 
in an arbitrary gravitational field is very different in that the non-Hermitian Hamiltonian is 
used in that work, in deep contrast to the explicitly Hermitian one (12.211) . 

Since the metric tensor is symmetric, it can be brought to the form diagonal in spa- 
tial coordinates. For example, the Kerr metric in the both spherical and Boyer-Lindquist 
coordinates belongs to this form. 

A spatially diagonal metric tensor can often be reduced to an isotropic form by an appro- 
priate transformation of spatial coordinates. In spatially isotropic coordinates, W a b = W5 a b 
and the final form of the line element is defined by 

ds 2 = V 2 c 2 dt 2 - W 2 5 ab (dx a - K a cdt) (dx b - K b cdt). (2.23) 

Evidently, none of the two transformations changes the temporal coordinate. Since only 
the spatial coordinates are transformed, the transformations do not change the physical 
frame (particle rest frame) used for a definition of the three-component physical spin. 

For isotropic metric (12. 23[) with the Schwinger gauge, the exact Hermitian Dirac Hamil- 
tonian reads j^] 

U = f3mc 2 V + °- [(a ■ p) T + JF{a ■ p)\ 

c rlc 
+ - (K ■ p + p ■ K) + — ( V x K) ■ S, (2.24) 

where 

7 = V/W. (2.25) 
It covers the general case of a spin-1/2 particle in an isotropic stationary metric. 



III. THE FOLDY-WOUTHUYSEN REPRESENTATION 



Let us now derive the FW Hamiltonian for stationary spatially isotropic metric f)2.23p . 
We perform the FW transformation of the exact Dirac Hamiltonian (I2.24p with the help of 
the method developed in Ref . [3] . Since the gravitational field is supposed to be strong, we 
do not make any approximations for the functions V, W, K a , and the only small parameter 
is the Planck constant h. In the FW Hamiltonian, we retain all the terms of the zero and 
first orders in h and the leading terms of order of fr 2 nonvanishing in the both nonrelativistic 
and weak field approximations. These terms describe the gravitational contact (Darwin) 

n 

interaction |6J. The computations are straightforward, and the final FW Hamiltonian is 
given by 

Hfw — 'H-fw + ^-fw ' (3-1) 

where 



H 



(i) 

FW 



kmc 4 



-, [S • ($ x p) - £ ■ (p x $) + ftV ■ $] 



2e' 2 + mc 2 {e',V}' 

he 2 ( l 1 

|-, [E • (0 x p) - S • (p x Q) + hV ■ Q\ I , (3.2) 



Wfw = 2 (K-P + P-K) + • (V x K) - ^ ^ 



/(2) C 



1 ,{^ 2 ,S-Q}|(3.3) 



mW 2 + -c 2 {^ 2 ,P 2 }, 



^ = v^ 2 ), $ = j 2 vy, 



(3.4) 



Q = p X V (K p + p -Jf)- V (if p + p • K) xp-px(px(Vx K)) - ((V x K) x p) x p. 

(3.5) 

The equivalent forms of the latter quantity 

Q = p x (p x (V x K)) + ((V x K) x p) x p + 2p x (p ■ V)K - 2<5 ab (V a K)p 6 x p, 
Q = px V(p • K") - V(K • p) x p + p x (p • V)K - 5 ab (V a K) Pb x p (3.6) 



may be also useful. Eqs. fl3.1IV - fl3.6p (which in the limiting cases agree with our previous 



results obtained in Refs. jg, 7, |9|) belong to the principal results of this paper. 

The dynamical equation for the spin is obtained from the commutator of the FW Hamil- 
tonian with the polarization operator II = /3S and is given by 

dU i 



dt h 



[H 



FW 



n 



x s + n {2) x n, 



(3.7) 



where f^ 1 -* is the operator of angular velocity of rotation of the spin in the static gravitational 
field, 

and the contribution from the nondiagonal part of the metric is equal to 

= ^VxK--{— 2 l — -AF\Q\X. (3.9) 

2 8 {2e> 2 + mc 2 {e',V} 1 V ; 

The two different matrices appear on the right-hand side of Eq. (13 .7p due to the fact 
that fr 1 ) contains the velocity operator rather than the momentum one. Since the velocity 
operator is proportional to an additional /3 factor and is equal to v = f3cp/e for free particles, 
the operator when expressed in terms of v, also acquires an additional /3 factor j^. 

The corresponding semiclassical formulas describing the motion of the average spin are 
then explicitly given by 

^ = f]xs = (n (1) + f] (2) )x S , (3.10) 

e'{e+mc z V) 2e 

n< 2 > = ~ V x K - F 2 Q, (3.12) 

2 4e'(e' + mc 2 V) 

where, in the semiclassical limit, 



e' = ^m 2 (^V 2 + c 2 p 2 7 2 , Q = 2p x V(p • K) + 2p x (p ■ V)K. (3.13) 
FW Hamiltonian (13 .ip can now be expressed in the simpler form in terms of tt^: 
H FW = Pe'+ C -{K.p + p.K) + hi.^ + ^:.^ 



, ___, TJ TTT .V-^K^^.V-g^. (3.14) 



,fr 2 mc 4 r i _ j „n 2 c 2 r i 

~4~ [2e> 2 + mc 2 {e',V} 
It is instructive to compare the classical and quantum Hamiltonians of a spinning partic 
In order to do this, one can start from the classical Hamiltonian of a spinless particle 



21|: 



c. 



/ m 2 2 ~ o6 \V2 Oa Oa Ob 

n class = { ^H^J + + 9A0, ~g ab = g ab - (3.15) 

Substituting the components of the general metric (12. ip . we recast this into 



In order to take into account the spin correctly, this Hamiltonian should be completed by 



the interaction term s ■ f2, along the lines of the general discussion of the Ref. 

Haass = \/m 2 ^V 2 + c 2 5 cd F\ F b d w a w b + cK ■ tt + g$ + s ■ Q. 



22|: 



(3.17) 



In the general case, f2 includes the both electromagnetic and gravitational contributions. 

For the case of the stationary spatially isotropic metric (12.231) . we obtain the resulting 
classical Hamiltonian (equal just to cpo) 



U dass = ^m 2 c 4 V 2 + c 2 T 2 p 2 + cK-p + s-fl, 



(3.18) 



where terms dependent on the electromagnetic fields are omitted. It is satisfactory to no- 
tice that the quantum Hamiltonians (13. ID . and (13. 14ft agree completely with the classical 
Hamiltonian (I3.18p . 

The similarity of the quantum and the classical Hamiltonians naturally leads to the 
similarity of the quantum and classical equations of motion which we are going to discuss 
now. 

IV. CLASSICAL DYNAMICS OF SPINNING PARTICLES 

The dynamics of spinning particles in strong gravitational fields is described by the gener- 
ally covariant Mathisson-Papapetrou 23l |24j theory. A different approach was developed by 
Khriplovich and Pomeransky 22j for the noncovariant 3-dimensional spin defined in the par- 
ticle rest frame. In general, the analysis of the Mathisson-Papapetrou equations is a difficult 
task and various approximation schemes were developed for their solution. By neglecting 
the second order spin effects, the Mathisson-Papapetrou system is reduced to 25] 

DU a .„ 1 



dr 
DS a 

dr 



fa 
J m 

0. 



2m 



(4.1) 
(4.2) 



The Mathisson force in the right-hand side of ( 14. 1ft depends on the curvature R^p* of 
spacetime. The physical spin s is defined in the rest frame of a particle. Taking into account 
that the 4- velocity is 

Z7< = 7(4 + ^4), (4.3) 
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we can recast the Mathisson-Papapetrou system into the 3-vector form 

dv „ v(v • £) „ 1 „ , A AS 

-r = £ - v 2 ' + v x B + -/ m , (4.4) 

dr c z 7 

— = nxs, n = -B + ^——^ 1 4.5 

dr 7 + 1 c l 



by introducing the fields £ and £? 22 1 



- U'T/ = £\ -f/T/ = e\ c B c . (4.6) 

This represents a clear analogy between the gravitational and electromagnetic fields that 
makes it possible to speak of the gravitoelectromagnetic type effects. One should not confuse 



the objects H4. 6f) with the usual gravitoelectromagnetic fields [26|, \2j\ that are defined in the 
weak-field approximation and that satisfy the Maxwell-like dynamical equations which are 
derived from Einstein's gravitational equation. In contrast, the fields (14. 6 p arise in a purely 
kinematic context and are defined for arbitrarily strong field configurations. Note also that 
they depend on the velocity of the particle, unlike the usual gravitoelectromagnetic fields 
l^ . I27I ] . Nevertheless, their application is quite helpful because the equations of motion of 
particles and their spins, written in terms of £ and B, look very similar to the corresponding 
equations of motion of charged spinning particles in the electrodynamics. 

Let us calculate these fields explicitly. We start from Eqs. (12. lip and (I2.12p and use the 
notation from Ref. Multiplying the connection coefficients ( 12 .!![) - ("12. 12[) by the 4- velocity 
( Eg) , we find: 

V* r m = y + 7 + + C^) (4.7) 
U l r M = yQ^-^W^d c V. (4.8) 

As a result, for the general metric ( 12. ip . the gravitoelectric and gravitomagnetic fields ( 14. 6p 
read 

£a = ^(cQ^f b -C 2 W b ad b V), (4.9) 



B a = ^ y- -E a - -T v a + e abc VC bc d v d j . (4.10) 

The physical spin precesses ( 14 .5p with the angular velocity f2 that can also be written 
explicitly in terms of the connection jo), Q] 

= e abc ip (~ry£ + r^V/ c 2 ) . (4.11) 
11 



As compared with the quantities f2^ and d,^ which describe the precession of the quantum 
spin using the coordinate time, fi„ contains an extra factor dt/dr = U = 7/V in view of a 
different parameterization using the proper time. 

Let us calculate the classical precession angular velocity in the Schwinger gauge explicitly. 
Substituting (14. 7\\ and (I4.8P into (14. lip , we obtain the exact classical formula for the angular 
velocity of the spin precession in an arbitrary gravitational field: 



fi s = 1 ( -Tv a 



V V 2 



+ 



7 abc 



7 + 1 



e abc w\d f ,y v . 



1 r abc Q 

7 + 1 ^ 



(4.12) 



The terms in the first line are linear in the 4-velocity of the particle, whereas the terms 
in the second line contain the even number of the velocity factors. For the special class of 
stationary geometries (12.231) . this general formula reduces to the simpler expression: 



V 



-Vv xV 



1 

TP 



7 



v x W 



V x K 



7 



2c( 7 + 1) 



(v x V(« ■ K) + v x (v ■ V)K) 



(4.13) 



Equivalently, 



1_ 
V 



1 1 
v x $ v x y 

JV(7 + 1) 2T 



V x K 



1 



2c( 7 



(v x V(w • K ) + v x (v ■ V)K) 



Since 



mc V\ 1 



W 2 m*c 



2r2 



mc 2 V'y, 



p = mW^v, 



(4.14) 



(4.15) 



we conclude that the classical equation of the spin motion (14. 5 1) agrees with the quantum 
equation ( 13. 7ft and with the semiclassical one (13. 10D . Thus, the classical and the quantum 
theories of the spin motion in gravity are in complete agreement. This is now verified for 
the arbitrary strong field configurations. 



V. EVOLUTION OF HELICITY IN AN ISOTROPIC METRIC 

Let us analyse the semiclassical evolution of the helicity of a particle propagating in 
a strong gravitation field. The helicity describes the spin orientation with respect to the 
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direction of particle's motion. As usual, one means the orientation of the 3-component spin 
defined in the particle rest frame (physical spin). The particle motion is defined by the 
trajectory that shows how the contravariant spatial world coordinates of the particle change 
in time. Evidently, the particle motion can be correctly characterized by the evolution of 
the contravariant world velocity or the unit vector in its direction. As a result, one can 
unambiguously define the helicity as a projection of the 3-component spin (pseudo) vector in 
the particle rest frame onto the direction of the unit vector along the contravariant velocity 
in the world frame. Thus, the helicity should be defined as 



where U = {U 1 ,U 2 ,U 3 } and V = U/U°. The investigation of the helicity is simplified 
when particle's trajectory is infinite. In this case, one can apply the fact that the vector 
N — U/U coincides with the vectors v/v and p/p at the initial and final parts of such 
a trajectory of the particle because of the very large distance to the field source. Here 
v = {y , v 2 , v 3 } is the velocity in the anholonomic coframe (14.31) and p = {— Pi, — P2, —P3} is 
the covariant momentum entering the classical and quantum Hamiltonians. For the isotropic 
metric under consideration, p a = e°pg + e b a pi = Wpa and p a /p = Wife/ (\W\v) = v$/v. As a 
result, the change of the helicity on the whole trajectory can be given by A£' = A£, where 
(' = (s/s) ■ n, n = v/v = p/p. Evidently, (' = cosx, where \ is an angle between the s 
and n vectors. 

In the general case, the evolution of the covariant momentum operator is defined by 



In the discussion of the semiclassical dynamics of the particle, we can neglect the small 
spin-dependent force which enters into the equation of motion with an additional factor h~. 
Thus, in the semiclassical approximation 



C = (a/ 8 ) ■ (U/U) = (s/s) ■ V/V, 



dp % 

~dT ~ h 



[H FW ,p\. 



(5.1) 



dp 1 1 dp 



dt p 3 dt ' 



and the dynamics of the unit vector n is given by 




(5.2) 



Therefore, the vector n rotates with the angular velocity 

- 1 dp 




(5.3) 
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and one can add a quantity nn with an arbitrary factor k. Using the Eqs. (I3.4p . (13.141) . 
dUSD, dEU), O, we derive 



c 2 f 

= ~ -py 2 n x ~ i^j; n x Q — cn x V(n ■ K). (5.4) 

The nabla operator does not act on n. 

We can add the quantity -n(n ■ (n x K)) to £c and use the identity 

2n x V(n • K) = n(n ■ (n x K)) — VxK + nx V(n ■ K) + n x (n • V)K. 

As a result, Eq. (15. 4p is recast into 

2 

u> = — n x $ n x C? + - VxK fn x V(n ■ Jf ) + n x (n • V)K1 . (5.5) 

The angular velocity of the spin rotation defined by Eq. (13 . 1 2[) can be expressed in the 
form similar to Eq. (14.141) : 

v v _ 
Q = r n x $ n x Q 

TV (7 + 1) 27 

+ - V x K - C<CI - ^ [n x V(n ■ K) + n x (n ■ V)K\ . (5.6) 
2 27 

Therefore, the vector of spin rotates with respect to the momentum direction, and the 
angular velocity of such rotation is 

c 2 c 

o = CI - uj = — n x <& + — [n x V(n ■ K) + n x (n ■ V)K] . (5.7) 

The same result can be obtained for an arbitrary metric by the purely classical method 
considered in the previous section. Indeed, let us use classical equations (14. 4p . (14. 51) and 
neglect the relatively small Mathisson force / m . Since 

^ = — - n x ( n x ^ ] , (5.8) 
dr v \ dr J 

Eq. (14 .4p leads to the following formula: 

dn , 

— = u'xn, 5.9 
dr 



where 



u'=-B+\vx£. (5.10) 
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Consequently, we find that the spin vector rotates with respect to the momentum direction 
and the angular velocity of this rotation is 

o=-(Q-u;') = -V^4^. (5.11) 



The gravitomagnetic field 3 does not influence particle's helicity [31] . This conclusion is 
valid for an arbitrary gravitational field when the effect of the Mathisson force on the particle 
motion is neglected. S is proportional to j/V. The use of the world time shows that any 
gravitational field does not affect the helicity of a particle with a negligible mass. 

The covariant Dirac equation can be used and the FW transformation can be performed 
for a massless particle. However, the interpretation of the results is a serious problem. The 
3-component spin is defined in the particle rest frame which cannot be introduced for such 
a particle. 

For stationary isotropic metric (12 .23j) . the gravitoelectric and the gravitomagnetic fields 
(1491) and (OUjl reduce to 

v~ [w daV + nr ~ + d{aKb) ~ ) ' (5 - 12) 



This demonstrates the perfect agreement between the quantum and classical analyses, cf. 
Eqs. ([57fD and ( IBTTll . 

As a final remark, we should mention that the results presented in this section cannot, 
generally speaking, be applied for the particle moving on a finite trajectory because the 
replacement of ( by (' is not always possible in this case. 



A. General noninertial frame 



A general noninertial frame is an important application of the results above. This frame 
is characterized by the acceleration a and the rotation uj of an observer. The exact metric 
of the flat spacetime seen by the accelerated and rotating observer has the form ( I2.23P and 
is given by Eq. (12. 1\ . 

We can apply our general results (13.71) ( |3T9l) . (13. 14ft for the metric ( 12 .2D of the flat space- 
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time. For the FW Hamiltonian we then find 



r ( a ■ r N 



^° = ? 1 1 + — > V m2(A + ° 2 P 2 - w • Z ' (5-14) 



2 IV c 2 

n« = QXp , o( 2 ) = 7 = v/m2c4 + cV . 

mc 2 (7 + 1) mc 2 

Here Z = r x p is the angular momentum operator. There are no any terms of order of 
h 2 nonvanishing in both the nonrelativistic and the weak field approximations. We also 
emphasize the absence of the spin-dependent Mathisson force. Let us stress that Eq. (I5.14p 
is derived for the strong kinematical effects when the ratios \a ■ r\/c 2 and \ui x r\/c are 
not small. Nevertheless, Hamiltonian ( I5.14p coincides in the case of uj = with the result 
6| obtained for the special case of the accelerated frame in the weak field approximation 
\a ■ r\/c 2 1. When a = 0, we reproduce the exact FW Hamiltonian deduced in Ref. 
for Dirac particles in the rotating frame. Eq. ( 15.141) and the corresponding Hamiltonian 

□ n 

obtained in Ref. [7j generalize the approximate expression derived by Hehl and Ni |10| . 
The equivalent form of "Ho reads 

nil \ iftP a P , /r 1 r\ 

fio = P \ mc +ma-r) r y u> ■ I. (5.15) 

2mc z 7 

The metric under consideration is non-Minkowskian at any parts of finite and infinite 
trajectories of particles. Therefore, we need to use the contravariant 4-velocity U" 1 or the 
world velocity V = cU /U°. Evidently, V/V = U/U. The world velocity operator is given 
by 



it h l • 2 |V c' I ' + p2 

As a result, the semiclassical approximation for the world acceleration operator 

w = IF = H [UFwy] 



reads 



. = -a(l + ^)-,x(2V + .xr) + (a ■ (2V + u x r)) + ^ x 
c 2 J c l + a ■ r 



The exact semiclassical Eq. (I5.17P agrees with the particular results obtained in Refs. 6| 
and |3] (which are approximate for the acceleration a). 



1(3 



The expression for the angular velocity of rotation of the vector N = V/V is similar to 
Eq. Q: 

-"•^"T (5 ' 18) 

Explicitly, 



iV 

wv = — 2oj + — x 



a • r\ . , fa - (2V + w x r)) , 
-a H — - uj x (w x r) + ^ — (u> x r) 



. (5.19) 



c z / c z + a ■ r 

The exact semiclassical expression for the angular velocity of spin rotation can be written 
in terms of V: 



a x (V + uj x r) 



UJ. 



(5.20) 



i + ^ + V(i + ^) -F^)' 

To consider the helicity evolution, it is sufficient to examine three specific cases. If 
\uj x r\ <^ c, \a ■ r\ c 2 , one may retain only the terms linear in uj, a: 

V x a 



Oy = fl — UJy = UJ 



7 V 2 



When a = (a purely rotating frame), we find 

V x (uj x (uj x r)) 



Oy = UJ 



V 2 



(5.21) 



(5.22) 



The presented solution for the rotating frame is exact. 

When uj = (an uniformly accelerated frame), we obtain 



Vxa 



1 



a ■ r 



Vxa 



1 + ^ 

cr 



2 _ y2 
73" 



(5.23) 



As follows from the Eq. f)5.16p . V/c = 1 + a • r/c for the ultrarelativistic particles with 
negligible masses. For such particles, we find Oy = 0. Therefore, their helicity remains 
unchanged in the uniformly accelerated frame. 

We can conclude that the motion in the general noninertial (arbitrarily rotating and 
accelerating) frame leads to the change of the helicity even for the ultrarelativistic particle 
with a negligible mass. A similar effect takes place in a gravitational field when a particle 
trajectory is finite. 

Another possible application of our general results is a rotating massive thin shell which 
metric obtained by Brill and Cohen 28| has the form (I2.3p .( l2~23|) . We will analyse this 
elsewhere. 
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VI. SQUARED DIRAC EQUATION AND THE EQUIVALENCE PRINCIPLE 



The electromagnetic and gravitational contributions to the covariant derivative (I2.10p 
manifest an obvious similarity of the electromagnetic and gravitational effects. In this section 
we further clarify this similarity by analysing the squared Dirac equation. 

The commutator of the covariant derivative (I2.10p reads 

D . D . _ DjD . = % A F .. + % -a aP R l3 ^. (6.1) 

Here Fij = diAj—djAi is the electromagnetic field tensor, and Rij a ^ is the Riemann curvature 
tensor. 

Acting with the conjugate Dirac operator {ihy a D a + mc) on (12. 9ft . we find the squared 
Dirac equation 

[ r Wg ii D i D i - lr- i ' ! - m 2 c 2 ) ip = 0. (6.2) 
Substituting Eq. (16. ip . we find 

(-ffgVDiDj - ^ 7 Vi ? « + J 7 [ V HolftRi^ ~ m 2 c 2 ^j j> = 0. (6.3) 



Expanding the product of gamma matrices, we derive the equation 



(-H 2 g ij DiDj - ^ a* Fa, + ^R- m 2 c 2 ) $ = 0, 



(6.4) 



where F a p = e^e^Fij are the tensor-like electromagnetic field coefficients. 

The special form of Eq. (16 .4p for the Dirac particle in a gravitational field has been 



obtained in Ref. 



29]. 



Explicit computation yields 



Since 



- h 2 g ij DiDj = TrVi - -a^ {tt\ T iaP ) + ^-a^a^r^V ilMV . (6.5) 

4 16 



we get finally 



ttV* - (^F Q/3 + m$ a/3 ) + - i? + — (2P a/3 r^ + ie^P^r^s) - m 2 c 2 



^ = 0, 
(6.6) 
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where 

= ^ T iaP ] , 75 = -iTVTY- (6-7) 

In the semiclassical approximation, n 1 = mil 1 , and $ ai g coincides with the spin (and 
momentum) transport matrix in a gravitational field (see Ref. js]) and with the tensor-like 
coefficients ^ a p\u x 22[. It is analogous to the electromagnetic field tensor and leads to 
the Dirac gyro-gravitomagnetic ratio g gra v = 2 in perfect agreement with the equivalence 
jrinciple which is also manifested in the interaction of spin with gravity j^o)]. This means 



3l| the absence of both the anomalous gravitomagnetic moment and the gravitoelectric 



dipole moment which are gravitational analogs of the anomalous magnetic moment and the 
electric dipole moment, respectively. 

Eq. (16 .6p explicitly shows a similarity of the Dirac particle interactions with electromag- 



netic and gravitational fields. This similarity is caused by the similarity of the motion of 
spinning particles in any external classical fields shown in Ref. 



VII. CONCLUSIONS 



In this paper we have studied the quantum and classical dynamics of spin 1/2 fermions 
in the strong fields. We derived the Hermitian Dirac Hamiltonian (I2.2ip that describes the 
fermion in arbitrary electromagnetic and gravitational fields in the Schwinger gauge. Apply- 
ing the FW transformation, we constructed the respective Hamiltonian ( 13. ip for an arbitrary 
isotropic metric. This allowed us to obtain the operator equation of spin precession (13. 7p 
and its semiclassical limit (I3.10p . As a specific application, we proved that the equations of 
motion in the uniformly accelerated rotating frame obtained earlier for the small accelera- 
tion are in fact valid for arbitrary large acceleration. We also derived the classical equation 
( I4.12p for the spin precession in the general gravitational field (in the Schwinger gauge) 
and observed its consistency with the semiclassical limit of quantum equation, in complete 
agreement with the equivalence principle. We confirmed this conclusion by deriving and 
analyzing the squared Dirac equation. The general evolution equations for the helicity were 
obtained in quantum and classical pictures. Their application to the dynamics in the general 
noninertial frame revealed the change of the helicity even for the ultrarelativistic particle 
with a negligible mass. 
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